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We present details of the phase diagrams of fermionic systems with random and frustrated interactions emphasizing 
, ^ ■ the important role of the chemical potential [i. The insulating fermionic Ising spin glass model is shown to reveal 
different entangled magnetic instabilities and phase transitions. We review tricritical phenomena related to the strong 
correspondence between charge and spin fluctuations and controlled by quantum statistics. A comparison with the 
diluted Sherrington-Kirkpatrick spin glass and with classical spin 1 models such as the Blume-Emery-Griffiths model 
is presented. We present a detailed analysis for the infinite range model showing that spin glass order must decay 
— i discontinuously as fj, exceeds a critical value, provided the temperature is below the tricritical one, and that the 
zero temperature transition is of classical type. Replica permutation symmetry breaking (RPSB) of the Parisi type 
describes the thermal spin glass transitions, along with modifications of the SK-models Almeida-Thouless-line. RPSB 
occurs in any case on the irreversible side of the (fermionic) Almeida Thouless lines and hence at least everywhere 
within a fermionic spin glass phase. Although the critical field theory of the quantum paramagnet to spin glass 
transition in metallic systems remains replica-symmetric at T = 0, with only small corrections at low T from RPSB, 
the phase diagram is affected at 0(T°) by RPSB. Generalizing our results for the fermionic Ising spin glass we consider 
aspects of models with additional spin and charge quantum-dynamics such as metallic spin glasses. 
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I. INTRODUCTION AND SUMMARY 



In this paper we concentrate on the magnetic phase diagram of the fermionic Ising spin glass and some of its related 
models, including a metallic extension. While in the preceding paper (denoted by I and the present one by II in 
the following) we showed the particular importance of Parisi replica symmetry breaking (RPSB)Era for the T = 0- 
and low-T properties of fermionic correlations we now consider mainly thermodynamic aspects. A fascinating and 
challenging selfconsistency problem shows up: knowledge of the phase diagrams is necessary to see where RPSB must 
be taken into account (no matter how complicated replica-symmetric calculations have already been) and finally 
the phase diagram, as paper I indicates, depends partially on RPSB. With respect to the absence of RPSB outside 
spin glass phases and in zero magnetic field, the quantum spin glass models of all colors are fortunately quite in the 
tradition of classical spin glassesHu. This helps to gain information for and from those outside areas without the 
complications of RPSB. 

In II, we present details of the tricritical behavior observed recently in thermal and in quantum-critical .transitions 
of both insulating and metallic model versions. A review of previous results for tricritical phenomena^ is joined 
with new details obtained for the discontinuous regime at higher chemical potentials, for low temperature solutions 
in general, and for instabilities inside the spin glass phase which are not primarily linked to replica permutation 
symmetry breaking. 

Considering fully frustrated magnetic interactions it is natural to expect spin glass order to play a major role in 
the phase diagram, and it will thus be surprising to find other types of transitions. Magnetic phenomena of these 
fermionic models, which are naturally described in the grand canonical ensemble, react strongly to the variation of the 
particle pressure controlled by the chemical potential. Charge and spin fluctuations are forced by quantum statistics 
to cooperate, which leads to the tricritical phenomena and to further instabilities in the paramagnetic regime as well as 
inside the spin glass phase. In addition to the conventional scenario of spin glass transitions (left aside the possibility 
of reentrant behavior) an unconventional critical line, linked to the creation of additional metastablc solutions, passes 
through the tricritical point and extends into the entire spin glass and paramagnetic phase. 

While the analysis of fermionic spin glasses deals only with real chemical potentials in first place, it is also clear 
that all so far considered quantum spin glasses (a prominent example being the transverse field Ising spin glass) can 
be viewed as fermionic spin glasses with properly chosen imaginary chemical potentials. This choice and the phase 
diagram of the model, comparing for example spin i- and spin 1-models, depends on the spin quantum number. 
Thus it becomes evident that understanding the phase diagram of fermionic spin glass models within the whole plane 
of complex chemical potentials will be important. We note that this point of view also applies to the twodimensional 
Ising model, recalling the success of other representation in terms of Majorana fermions or what has been called Ising 
fermionsErLl. The application to experiments which are described by fermionic systems with random interactioa-js very 
likely a wide open fielaI3. Yet a large number of experimental results already exists and we cite here a fewEm. The 
natural field of application includes HighT c superconductors, heavy fermion systems, semimagnetic semiconductors 
etcetera. 

One must be concerned in general with low temperature behavior near (T = 0) quantum phase transitions, which 
is very susceptible to quantum dynamical effects, but also interesting types of (asymptotically classical) thermal 
tricritical transitions which mix span- and charge-fluctuations occur. Interference of spin glass features in charge 
response and transport propertieai-3 is of central interest. 

The symmetry classification of the concerned QPT's turned out to be very different from that of thermal phase 
transitions and also does not appear to resemble, for example, the T = QPT universality classifications known for 
Anderson localization. i_ . ._. 
In early papers, using the spin-static approximation, the fundamental question had been raisedllTO whether tunneling 
through energy barriers in contrast to thermal hopping over them might distinguish quantum spin glasses from classical 
ones. The quantum-mechanical image of this type of built-up of zero temperature spin glass order could have been 
given by a Parisi order parameter function q(£, x) in place of the known classical one q(T, x) with £ being a parameter 
driving the system at T = through the transition. In the preceding paper, cited henceforth as paper I, we have shown 
that this is not the case. Instead the quantum dynamical image of Parisi replica symmetry breaking is extremely 
important even at T — but in a totally different way. It determines elementary features of the fermionic correlations. 
The strength of these effects as observed at half-filling and the qualitative changes entrained by them leave open the 
question in which way the magnetic phase diagram is affected by Parisi RPSB. This is very hard to answer for the 
moment, since the analysis may require a sufficiently good approximation of the Parisi function at T = (a fact that 
has not been necessary at half-filling) . On the other hand our analysis of paper I did not show any sign of a limitation 
of RPSB- effects to half-filling. On the contrary unsurmountable difficulties at finite order of RPSB in deriving stable 
homogeneous saddle-point solutions away from half-filling could possibly be cured by use of the full fermionic Parisi 
solution for low temperatures. This requires however the knowledge not only of the Parisi function q(x) at low Td'tl 
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but also of quantum-dynamic fermionic Green's functions as described in the preceding paper. 

Quantum spin glasses form indeed a link between the general statistical theory of spin glasses and randomly interacting 
systems on one hand and of interacting many fermion systems on the other. Beyond this secondary role they also 
assume their unique place. There are clear-cut differences between for example T — spin glass transitions and 
other quantum phase transitions, be they magnetic or electronic like the Anderson localization. As for their thermal 
transitions they can, despite some relationships, as well be clearly distinguished from random field systems for example. 
Charge correlations and fluctuations in metallic spin glasses are of particular interest, since electronic transport uses 
the part of the Hilbert space which is spanned by nonmagnetic states. Magnetic transitions alter nonanalytically the 
occupation of these states which in turn leads to nonanalytic charge fluctuationsO. These secondary critical phenomena 
comprise for example the effect of non-Fermi liquid behavior ia the vicinity of a quantum spin glass transition^, for 
metallic quantum spin glass transitions had beeHr-given beforeES, emphasizing similarities and differences with respect 
to the one for transverse field Ising spin glasseaia. Charge variables are now included to describe the corresponding 
fluctuations. 
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II. TRICRITICAL PHASE DIAGRAM OF AN ISING SPIN GLASS WITH CHARGE FLUCTUATIONS 



One crucial feature of fermionic spin glasses is the intimate connection of spin and charge degrees of freedom. Quantum 
statistics with the presence of nonmagnetic states (single and double occupancy) opens the possibility of charge 
fluctuations which thermally dilute the spin system. The simplest model displaying this type of behavior is the 
SK-model on a fermionic space with four states per site, denoted by ISGf and defined by the Hamiltonian 



(1) 



with spins <x; = a Q' i.g ^fi r , particle number operator m = 9} a ^i, a and Gaussian distributed exchange integrals 
Jjj with variance J 2 H3. The fermionic field operators obey the usual commutator relations {^>i a , ^jp} — and 
{^la> ^jp} ~ As a guide to the global phase diagram we study an exactly solvable infinite range version of the 

model but the formulae obtained in this subsection may equally well be considered as the saddle point approximation 
for an interaction with finite range. We recall the free energy (thermodynamic potential) for arbitrary order K of 
Parisi replica symmetry breaking 2 ^ 
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cosh{/3H) + cosh(/3/z) exp(--/3 2 J 2 (g - qi )) 



(2) 



The abbreviating notation with upper index G denote normalized Gaussian integrals over z, ie dz eXP ^ 2 ^ ' ; 
and q(x) coincides with the Parisi order parameter function (the replica diagonal part q is not contained in it). The 
fields z v are needed for decoupling the m u x m„ blocks of the Parisi matrix Q and the effective field H is the sum 
of the external field h and the decoupling fields JY^=i \f*b> ~ Qv+i z v Notice the appearance of the Parisi variables 
q u and the additional q which lies at the heart of the following discussion. The Edwards - Anderson order parameter 
qEA — hm^oo < Si(t)Si(0) > describes the freezing of spins in the spin glass phase and is given by the Plateau height 
q(l) of the order parameter function, whereas the replica diagonal q = [< <j a a a >] av is related to the average filling 
factor [v] av = [< ^Jj^q >]av by [v] av = 1 + tanh(/?/i)(l — q). The last relation is exact even in the case of replica 
symmetry breaking. 



A. Spin glass transition and unusual critical line 



To understand the effects of charge fluctuations and to gain insight into the global phase diagram it is sufficient to 
consider a replica symmetric approximation of eq.(||) with q{x) = q constant 



/ = l/4/3[(l - qf - (1 - qf] - Tin 2 - (J, — T f hi{C M (z)}, 

J Z 

C^z) = cosh[/3 J ff(z)] + cosh(/fy) exp[-l/2/3 2 (g - q)} 



(3) 
(4) 



where H(z) stands for J^fqz + h. Differentiation with respect to the saddle point variables q and q yields the 
corresponding selfconsistency equations 



smh 2 [(3H(z)} 
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q = 
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whereas differentiation with respect to h and \i yields the equations for the magnetization and the average filling [v] av 
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sinh(!3H{z)) 

CM (n 

[u] av = [< ^ a >} av (8) 

- l + J z C,{z) (9) 

Phase transitions are signalized by vanishing masses of the order parameter propagators which in the saddle point 
formalism are given by second derivatives of the free energy. On the other hand a positive mass for q and a negative 
one for q guarantees stability. In the paramagnetic regime the relevant expressions are 



^ = -l'^' 2 7(l -<li -0 (ID) 



= ^[-i + /?Y]<o (11) 

This approach correctly tests critical fluctuations in the paramagnetic regime, in the following sections a more complete 
analysis including replica symmetry breaking fluctuations will be carried out using the a»«roach of Almeida and 
ThoulessEH. A similar system of coupled stability conditions .was found for the BEG - mode£3 and for a SK - model 
with crystal fielcO, the case of half filling was considered inc3. Analyzing ( pd| ) as an equation we get the condition 
T c = q c for the transition to spin glass order and solve (||) for 

/i cl (T) = Tcosh- 1 [(l/T- l)exp[l/(2T)]] (12) 



In the same way we calculate the curve corresponding to critical charge density fluctuations 

^c 2 (T) = Tcosh- 1 ^ 1 T y/1 ~ 8T ^ 2 expf 2 11 (13) 

^ y ' L 8T 2 PL lT\/l-8T 2 

The graphs of these two functions are displayed in figure 1; they cross each other at /x = .958, T = .35 and have 
a common tangent point at fi = l/3cosh [2exp(3/2)] = 0.961056, T = 1/3. There seems to be a region where 
the diagonal elements of the Q - matrix (e.g. charge fluctuations) become critical before the off - diagonal elements 
which may indicate some type of phase transition. To shed some light upon this problem it is instructive to solve 
the selfconsistency equation (^) for q = 0. For T — ► oo there clearly exists only the physically meaningful solution 
q = 1/2. The solution stays unique until one crosses (coming from large temperatures) the [i c i line the first time. 
At this temperature two new solutions show up - one maximum and one minimum of the free energy. However, the 
physical solution still has a lower free energy than the new solutions and in addition to that the new solutions violate 
the stability condition (|Tl|). For .95 < /i < .962 one crosses the fi C 2 line a second time when lowering the temperature. 
Now something more severe happens: the maximum which at the first crossing coincided with the second minimum 
moves to the physical minimum, merges with it and then disappears. This doesn't have any consequences as the 
system is in the ordered phase at this temperature anyway and playing with paramagnetic solutions seems to be of 
academic interest only but there is one important exception: the point (/i = .962, T = 1/3). At this point the physical 
(paramagnetic) solution is unstable to both diagonal and offdiagonal fluctuations, a situation that deserves further 
investigation. The second surprise is that the fi c i(T) curve which at least for small fj, corresponds to a second order 
phase transition to a spin glass phase turns around at fi — .961056 and thus leaves us with the question about the 
nature of the low temperature state of the system. To gain some insight into this question we look at exact low 
temperature solutions to the selfconsistency equations (||) and (j|). 
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B. Exact low temperature solutions 



To analyze solutions of the saddle point equations at low temperatures it is useful to introduce the abbreviation 
x = cosh(/3/x) exp[— l/20 2 {q — q)] and to integrate the difference of eqs. (||) and (^|) 

i —— i i j.-^ i i i . / . i i .1.1. 

(14) 



(1 



1 + cosh(/3 v /gz)a; 
(cosh(/3/i) + x) 2 



T 



/2nq 
2T 
fhtq 



1 + xcosHy) = 
, (x + cosh(y)) 2 v ' 



+ 0{T 3 ) , 



(15) 



surprisingly enough the result is independent of x. 

The leading low temperature contributions to the q- equation can be integrated exactly. Its relevant asymptotic 
behavior, needed for the following discussion, is extracted by 



xT 

V^M J-c 
2xT 



dy 



1 



cosh(y) 



(16) 



f for x 



/2nq \ \ ln(2x) for x — > oo 
The behavior of x in the low T limit is determined by the sign of — P/2(q — q) 



IMI 



For < 



vanishes and we obtain q = 1 — ^/?Texp[/3(|/x| 



7m^> q 



1 - J^T. For a 

V 7T r 



/2tt 



x stays finite when T goes to 



zero and one still has q(T — 0) — q(T = 0) = 1. For still larger values of the chemical potential this value decreases, 
since now x grows exponentially and eqrtlq) turns into 



2xT 1 



2nq x 



ln(x) 



i 



/2irq' 



(17) 



The last equation is of fourth order in ^fq and can have several solutions. In the limit [i — > oo however the solution 
becomes unique and an expansion in inverse powers of the chemical potential seems useful. One finds 



9 = ^(1/^ + 1/^ + 0(1/^)) 



(18) 



We now compare the free energy of the paramagnetic solution with that of the ordered one. 

Below the /i C 2-line there exist three paramagnetic solutions, two minima and one maximum of the free energy. The 
maximum with the lower free energy (q — 1) is unstable with respect to spin glass order, hence the physical solution 
is 

q para = 2e-^ + 0(f3 2 e~ 2 ^)- (19) 
The corresponding free energy at zero temperature is fdis — — whereas the free energy of the spin glass solution is 

1 

4 A 

x 



f ordered = -f3[-2(q - q) + q 2 - q 2 ] 



/' 



(20) 
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cosh(/3 yfqz) 
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ln(cosh((3 y/qz)] 
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\/2vrg Jo 
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d.v 
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dz z exp( — — ) 
1 



1 



"(1/2 + 
A* 

1/1 1 X 



x + cosh(y) 
2~ 



1/4) 
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Hence at least for large values of the chemical potential this replica symmetric analysis indicates a paramagnetic 
ground state. On the other hand we knownJ that for half filling the effect of the fermionic space is to lower the 
temperature of the continuous spin glass transition by a factor of 0.6767. Hence a transition from spin glass to 
paramagnet must take place on the zero temperature axis. In order to locate it we solve Eq.(|l7j) for all values of the 
chemical potential and compare the free energy of the different solutions. Spin glass order is realized up to a chemical 
potential /i — .900, at this value we find a first order transition with a jump of the order parameter from q = .603 to 
zero. By means of the relation [v] av = 2 — q (valid in this form only at T = 0) one finds phase separation for fillings 
1.34<i/<2. nn 
However, a stability analysis of this replica symmetric solution according to the scheme of de Almeida and ThoulessEi'c£l 
reveals two negative eigenvalues of the Hessian matrix instead gf the expected one negative value indicating just the 
instability towards replica symmetry breaking. As described intj for spin glass problems, one has to choose the stable 
solution with the lowest free energy, hence the system undergoes a first order transition from half filling to the com- 
pletely filled paramagnetic state at /i = -j^J ■ The analysis of RPSB in I indicates that this region of incompressibility 
becomes smaller with increasing order of symmetry breaking. Therefor one may expect that for infinite order RPSB 
the fermion filling increases continuously from one with increasing chemical potential until it makes a finite jump at 
a special value of fi. 
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FIG. 2. Qualitative phase diagram of the ISGf in (/i,T)-space calculated in the replica symmetric approximation 

Now we have gathered important information on the global phase diagram of the infinite-range fermionic spin glass 
model. At /i = (half filling) a continuous transition takes place, for increasing fj, the transition temperature is 
lowered. At zero temperature the transition is first order in both the spin and charge density, therefore it is natural 
to look for a tricritical point at which the transition changes its character from 2nd order to 1st order. The point 
T = J/3, /j, = 0.961056 is a good candidate for tricritical behavior since charge- and spin- fluctuations become 
simultaneously critical there. 



C. Tricritical Behavior 



Resorting to different techniques we performed a detailed analysis of tricritical behavior in the ISGf. Mean field type 
quantities are most easily calculated from an expansion of the saddle point free energy around the tricritical point 
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given by 



= J/3 and ^ c3 



^arcosh(2 exp( — )) 



.9611J 



(21) 



This expansion reads 



TCP 



= M~Mc3 - ~2j( S ^~ J l( x ) dx ) 



J{(lr g g-r T 6T')6q + lsT[(Sq) 2 



dxq 2 (x)} --[ dx[xq 3 (x) + 3q(x) / dyq 2 (y)} 
o * Jo Jo 



(22) 



where Sq = q — qTCP,gJ = I 1 ~ ^c3 + 3(C 1 J — fi c s)ST as nonordering fielcc3, and <5T = T — T c ^. The constants 
are given by r g = ^-,Tt = 12(1 — fC -2 ) with £ = tanh(fi C 3/T C 3) ~ 0.9938, and /i C 3 as the characteristic chemical 
potential locating the TCP. The average filling factor corresponding to /^ C 3 is evaluated as [z^ou — 1-6625. Moreover, 
we derived a fluctuation Lagrangian for the tricritical and finite range ISGf, a Lagrangian of the same structure is 
obtained for generalized models (eg with a transport mechanism) at finite temperature by integrating out dynamical 
degrees of freedom 



L = - I d d x 



rK\ 



°-h 2 Y,Q ab 

® aa + Q aa (-^ 2 + u)Q a 



1 1 ' 

-rr'(vg Qb ) 2 - - Q aa Q bb - ^]T(Q aQ ) 3 



«3 



7VQ 3 - k 2 Q aa Q ab Q ha + fE (0 afc ) 4 ], 



(23) 



Here 4(^-) 



(0) 



(M)(0) 



^- and u' ^ 



denote the bare coefficients at tricriticality. One fourth order 
term relevant for replica symmetry breaking is kept. Replicas under or Tr' are distinct. The Q aa Q bb -coupling is 
renormalization group generated as in the metallic quantum spin glass, its effects will be discussed in a subsequent 
section. The simultaneous appearance of critical diagonal and off-diagonal field components is very unusual for 
classical thermal spin glass transitions and was so far known only to occur in special limitscl 

The tricritical behavior can be extracted from a quadratic approximation of the self - consistency equations. Since 
replica symmetry breaking is only generated by the y± - term in the free energy in the context of a quadratic 
approximation it is sufficient to consider the plateau height q{l) = qEA = <Z+ subleading corrections. 



From 



we get for q=0 



= gr g - r T ST 2 + 65TSq 
= 6q{5q -ST -q) 



-Sq 2 + 3q 2 



(24) 
(25) 



1 



5q dis =A{6T±\8T\ 
= A(ST±\ST\W) 



r g9 
VIST 2 



r JL" 
12' 



(26) 



Only the smaller solution (eg that with the - sign) corresponds to a minimum of the free energy; on the line Sfi = 
— 3(1/C — /z C 3)5T = — 0.1354(5T which is tangent to the critical curves fx c i(T) and fi C 2(T) g vanishes and Scjdis is 
proportional to ST. The statement holds true for a region close to the tangent where g ~ ST 2 as well, there the 
radicand is positive for g < — 12 ~ rT ST 2 . This result reproduces an expansion of the exact relation /j, C 2(T) around 



T = 







However, if g is of order ST or larger the solution becomes to leading order Sq = «/ ^ and thus displays a nonanalytical 
dependence on temperature and / or chemical potential. This type of crossover can also be seen in the free energy 

/*. - l\5T\"Wgn{8T) - W)(^ - r T ) (27) 
-A8sgn(ST)(sgn(ST) - W) 2 - l6(sgn(ST) - Wf] 
= \ST\ 2 - a G(^) 

The scaling form allows for the identification of the specific heat exponent a = — 1 and the crossover exponent <f> — 2. 
The crossover function Q (x) is regular for small x and has the asymptotic form 

~ 2 — a 

G(x) x * Qoo + (regular corrections) . (28) 

x — ► oo 

In the tricritical region the leading singularity in the free energy is given by 

4 3 

frcp^ng = —=(r g g) 2 (29) 

For S/i = we have g ~ ST and can read off the tricritical specific heat exponent = \ from above. 

The next step is the search for ordered solutions of the selfconsistency equations. Using q = Sq — ST one solves readily 

for 



Sqordered = ±^\5T\J -^j^ + r T - 3 (30) 



The solution exists for g < — - — (tt — 3) (see curve d in figure 2). The corresponding replica overlap q is positive 
for q = \5T\(—sgn{8T) ± | y — ^jf- + r T — 3) > 0. For positive ST the "+"- sign needs to be chosen, the condition 
q > is then equivalent to g < — (2i — r^)^— , which reproduces an expansion of /i c i(T) (see curve a in figure 2). 



Considering now the case ST < 0, we get the condition ±§y — Jflr +rr — 3 > — 1. The "+" - solution is always o.k., 

the "-" - solution only in the region with — — rx)^— < g < ~^—{tt — 3) (in the region between curve a and d). 
The saddle point energy for the physically meaningful " +" - solution is 



_ _2,_r g g _ r T ~3 ^ 

Jordered ^ j^-, 2 r ' \ 01 ) 



The energy difference f ordered — fdis becomes negative for g < — 7.9488(5T 2 (thermodynamic transition, shown in curve 
c in figure 2). In the tricritical region the solution of the selfconsistency equations is 



q = Sq = ~^-r g g (32) 

which yields the tricritical order parameter exponent fiz = \ and suggests that q and Sq act as order parame- 
ters simultaneously. From the fluctuation Lagrangian one reads off mass squared proportional to ST 2 and hence 
73 = 03 = a 3 = |. 



D. Replica Symmetry Breaking and tricritical Almeida Thouless - Line 



In contrast to crystal-field split spin glassesE3 the quartic coefficient j/4 of our free energy, Eq. (|22j) , is nonzero and one 
obtains the Parisi solution 

1 \ f for 0<x<x 1 

q(x) = < z v\ (33) 
(_ q\\) for X\ < x < 1 

The plateau height is found to satisfy 
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q(l) = Sq + 0(Sq 2 ). 



(34) 



Consequently, plateau and breakpoint scale like yjr] + 0(t) at the TCP, while linear r-dependence is reserved to 

T c > T C 3. Adapting the notation oiS we express our result for the irreversible response g(l) — f} q(x) ~ |r|^ A in terms 
of the exponent /3a3 = 1 for T — > T C 3 and /3a = 2 for T — ► T c > T C 3. For the Almeida-Thouless line at tricriticality 
we find 

H 2 80 2 /i c3 ,/3^c3 sss3/2 3/2 ri , 2 \ /or\ 

J 2 " = 8T r ~ — tanh (—))) t at + °( t at) (35) 

with = Tc3 ~Tat(.h) _ jjence we obtain the critical exponent $3 = | near T c3 , while = | for all T c > T c3 . These 
values do not satisfy the scaling relation 6*3 = with /3a 3 = 1 + (73 — cl§)/2. Along the lines described i this 
problem of mean-field exponents will be resolved below by the renormalization group analysis of the coupling j/4 of 
the finite-range and finite-dimensional ISGf. 



E. Stability of the solution in the tricritical region 

A first test for stability of the saddle point solutions is looking at the (matrix) of second derivatives of the free energy 
with respect to q and Sq. 



d 2 f 9 
dq 2 4 f 



9ST - -Sq + 0{q 2 ,6q 2 ) (36) 



^1 = 9(-5T + Sq - 2q) (87) 



In the tricritical regime 5q and q are leading compared to the temperature deviation. The stability condition < 

d 2 f • 

is satisfied both in the ordered and in the disordered regime, whereas is positive in the paramagnetic phase 
and negative (e.g. unstable) in the ordered phase. One may wonder whether the coupling between diagonal and 
offdiagonal fluctuations cures this problem or whether there is a new type of instability in addition to the well known 
AT - instability. . The SQ aa - fluctuations are vectors in replica space and play the role of the 5m a -fluctuations in 
the notation of Almeida and Thouless, the matrix entries A, B, C and D of AT are replaced by (to order ^~g) 

X = i/3 2 [l - l -(3 2 {< K) 4 > - < (O 2 > 2 )] 
27 

= ~Sq (38) 

Y = -/3 4 {< (a a ) 2 (a b ) 2 > - < {a a ) 2 >< (a b ) 2 >) 
= 

V = -P\< {a a ) 3 a b > - < (a a ) 2 >< a a a b >) 
= -54g 

W = -/3 4 (< {a a ) 2 a b a c > - < {a a f >< o b >< a c >) 
= 

P, Q and R control as usually offdiagonal fluctuations and are given by 

P = /3 2 [l- (3 2 {< {a a f{a b f >-< a a a b > 2 )] (39) 
= -M6q 

Q = -(3\< {a a ) 2 a b a c >-< a a a b >< a a a c >) 
= -27q 

R = -/? 4 (< a a a b a c a d > - < a a a b >< a c a d >= 

In the limes n — > 1 AT obtained five eigenvalues, one of them (P — 2Q + R) is not related to diagonal fluctuations, 
the other four merge in the replica limit to 
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A± = -(X -Y + P-AQ + 3R 

±y/(X -Y -P + AQ- 3i?) 2 - 8(V - W) 2 
= ^(X + P - AQ ± y/(X -P + AQ) 2 - 8V 2 

= ^(7± V2 9 - 3 4 ) (40) 
8 

The real part of these eigenvalues is positive and guarantees stability, the logarithm of the partition function stays 
real in spite of the imaginary part of the eigenvalues, it contains a factor 



v/(A + A_)(A+A_)(n-l) 

which is real in the replica limit. 



F. Tricritical behavior in the Spin 1 Ising model in transverse fields 



The quantum phase transitions at T — of metallic spin glasses and of the transverse field Ising spin | glass were 
recently described in a unifying quantum field theoretical way. The essential difference between the two being the 
marginal relevance of the quantum dynamical interaction of transverse field models in contrast to their (dangerous) 
irrelevance in the metallic case. Similar properties of the phase diagrams of spin 1 and of the ISGf suggest a com- 
parison between quantum extensions of these two models. The Spin 1 Ising spin glass in a transverse field will see its 
thermal tricritical point descend towards zero temperature as the transverse field approaches a characteristic value. 
Again, as in the metallic case, we approximate this value by a Q-static approximation, improving this result finally 
by a generalized Miller-Huse methodc3. 



G. Renormalization Group Analysis 



1. Tricritical Ising Spin Glass 



We performed a 1-loop RG calculation for tricritical fluctuations. At each RG step the mass of charge fluctuations 
5Q aa was shifted away. Introducing the anomalous dimensions fj and rj for diagonal and offdiagonal fluctuations one 
finds at one loop level the following RG relations (e = 8 — d) 

^ = (| - 11k£ + 16/si/sa + 6«§)r - «| (42) 



2(1 - k\)u - An\ + Ak x k 2 , 



du 
~dl 

dni e o dn 2 ,£ o 2 r. \ 

— = -m + 9k x — = (- + 6k 2 - k x + 16kik 2 )k 2 



-^- = (2+9^3, 

Above d = 8 the RG flows towards the Gaussian fixed point with mean field exponents, for d < 8 there is no 
perturbatively accessible fixed point for real k's. However, a preliminary analysis of the resulting strong coupling 
problem shows that there exists a solution with positive fj in contrast to the negative anomalous dimension typical 
of cubic field theories with imaginary coupling. The RG for the DIC 2/4 showed that its long-distance behavior is 
dominated by a K 4 -contribution (like inEj but) for d^ = 8 < d < 10. This leads to the modified MF exponent 

$3 = &/(d— 4), which satisfies the scaling relation #3 =J2Jj3/± 3 in d 1 ^ — 8 and reduces to the MF-result in 10 
dimensions. The dimensional shift by 2 in comparison withcJ is due to coupling t. 
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2. Metallic Quantum Spin Glass 



In the theory of quantum phase transitions time dependent fluctuations are treated on an equal footing with spatial 
fluctuationsEj. While for the Lagrangian (|2^) it was sufficient to consider only the uj — - component of the Q 

- fields, in the quantum case all low energy fields must be kept as they are coupled via the quantum mechanical 
interaction u J d,T(Q aa ) 2 . However, the value z = 4 of the dynamical critical exponent renders the u - coupling 
dangerously irrelevant and allows for a perturbative mapping of the critical theory on a classical problem, the Pseudo 

- Yang - Lee edge singularity. This field theory has only one cubic coupling which corresponds to K\ in eq(|2^). The 
comparison of the metallic quantum case with the thermal tricritical theory allows one to understand the nature of 
the strong coupling RG fixed point: whereas the spin fluctuations in the thermal second order regime are governed 
by a perturbatively accessible fixed point, the TCP and the quantum case are characterized by the combination of 
charge and spin fluctuations and a corresponding strong coupling problem. 

III. THE PARAMAGNETIC PHASE: MULTIPLE LOW TEMPERATURE SADDLE-POINT SOLUTIONS 

Complications due to variation of the chemical potential can also be seen outside the spin glass phase of the fermionic 
Ising spin glass. We have discussed (see also Fig.(|l|)) the unusual critical line ^(T), which crosses the spin glass 
regime and, passing through the thermal tricritical point, stretches into the regime up to arbitrarily large [i. It was 
seen that in the regime of large particle pressure the spin density becomes small and spin glass order is apparently 
absent. Despite the presence of the /^(T)- or, if one wish, T C 2(/i)-line, the absence of magnetic phase transitions 
caused by the totally frustrated magnetic interaction is predicted by the replica symmetric ordered solution eq. (ft7|) . 
The multitude of paramagnetic solutions below the T C 2 (/i)-line and a possible metastability in finite range models is 
discussed in this section. We take advantage of the fact that replica permutation symmetry will play no role. We 
consider however possible that spontaneous vector replica symmetry breakingcj will be involved. In the paramagnetic 
regime one finds three solutions for q, and hence for the susceptibility x = PQ too, at temperatures below T&ijJ-). 
This temperature is given by the inversion of Eq.([i3|) and is also shown in Fig. (5). 
At low temperatures these solutions are well approximated by 

q\ T<<T = | 2/X//3 + 0(l3~ 2 lnP) (43) 

These saddle-point solutions for q and the corresponding ones for the linear susceptibility are displayed for a few 
characteristic values of the chemical potential fi, starting from the tricritical value /z C 3 = 0.96125 J. The magnetically 
saturated solution (q — I) which would imply a Curie susceptibility \ ~ (3 is unstable with respect to spin glass order 
for all temperature. The high temperature solution, which decays exponentially at low temperatures and leads to a 
vanishing susceptibility in the zero temperature limit is the physical paramagnetic solution in the whole temperature 
range. It is not clear however whether a replica symmetry broken stable solution with a possibly lower free energy 
exists due to the presence of exceptionally strong bonds. The linearly decaying solution is a maximum of the free 
energy at all temperatures. For large values of the chemical potential the physical solution q = 2e~^ M has a much 
larger free energy than the Curie type solution q = 1. The two free energy minima are separated only by a very 
small barrier, hence we expect fluctuations to play a major role in finite range models. The instability of the lower 
minimum with respect to spin glass order makes it conceivable that due to Griffiths effects a replica symmetry broken 
ordered solution persists in the region of very low T for all values of the chemical potential. 



13 





(c) f3 (d) f3 

FIG. 3. Spin autocorrelation function shown for various chemical potentials /i — .96125J(a), J(b), 2J(c), and 10J(d) with 
three low temperature solutions of the paramagnetic saddle-point equation: the lowest one (indistinguishable over a wide range 
from the /3-axis on the given scale) is the only stable minimum of the free energy. 
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(H = Hcs = 96125J) 



(3 




(H = 33) [3 (i 

FIG. 4. Linear susceptibility corresponding to the three solutions of the paramagnetic saddle-point equation, given for \i=\, 
fi = 3, which are joined in the fourth figure (bottom right); the first one (top left) shows the solution near the tricritical chemical 
potential for comparison. The (sole) high temperature solution remains realized at all temperatures. 
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FIG. 5. Paramagnetic free energy as a function of q for n — 3 J and /3 = 10/ J (curve a),/3 = 20/J (b), and /3 = 50/J (c). 

IV. RELATIONS WITH THE BLUME-EMERY-GRIFFITHS MODEL FOR HE 3 - HE 4 MIXTURES; 

UNIVERSAL MEAN FIELD TRICRITICAL POINT 

For classical spin systems there exists a variety of different ways to describe diluted systems. Spin 1 models, where 
Si = may be considered to describe an empty site i, provide one example. A beautiful case was formulated by 
Syozi, who extended Onsager's solution to the twodimensional diluted Ising ferromagnet. Another example is given 
by the Blume- Emery-Griffith model intended to describe He 3 — He 4 mixtures, where the Si = state corresponds.-to 
a iIe 3 -atom on site i. Disorder effects added to the original BEG model have raised considerable interest recently^. 
But it is in fact the clean BEG model that shows the perhaps accidental but in any case most eyecatching resemblance 
with part of the fermionic Ising spin glass tricritical behavior. There is a surprisingly identical tricritical temperature 
T C 3 = \J of the clean BEG without S 12 ^-interaction (ie K=0) on one hand and of the ISGf on the other (the meaning 
of J being n.n. for the BEG and standard deviation of Jij in .case of the spin glass). Another surprising fact is 
however that it is the K ^ BEG-equation for < S 2 > (= 1 — x ir£3) which can be mapped onto the ISG/-equation 
for < (S a ) 2 >~ q in the disordered phase. Reconsidering the g-equation for H = m = q = from the selfconsistent 
i^Gy-cquations given above 



Up to now the tricritical phenomenon and the hrst order transition from paramagnet to spin glass have mainly 
been discussed for given chemical potential. In this picture a thermodynamical first order transition is observed for 
6/j, = — 0.13555T— 7.9488<5T 2 . On this line the free energies of paramagnetic and ordered solution are the same but the 
two solutions correspond to different fermion fillings of the system. This can be seen from the following consideration: 
On the first order line we have Sqdis = — 6.2649|<5T| and 5q rdered — 1.5133|<5T|. Via the relation 




(44) 



and comparing it with the BEG-equation for m — which reads 



1 - x = 1/[1+ -exp(0A)exp(-PK(l - x))] 



(45) 



one finds an identical structure in the variables q and 1 — x. 



V. PHASE SEPARATION 



1G 



5u = -Sq tanh(/3 c3 /i c3 ) - ST-^ (46) 

the phase diagram in (u,T) space in the neighborhood of the TCP is defined. 8v\ — 6.2522|(5T| limits the ordered 
region, 8V2 — — 1 .4779 |<5T| the disordered region. For intermediate fillings 8v\ < 8v < 8v2 both phases coexist. 



VI. THE FERMIONIC ISING CHAIN 



The remarkable similarities between the phase diagrams of the BEG-model without disorder and the one of the 
fermionic Ising spin glass can be taken as indicative for the fact that rather dilution than disorder is the source of 
the tricritical crossover from continuous to discontinuous phase transitions. Moreover the RG theory showed that the 
tricritical structure is also not very sensitive to spatial dimension. These impressions and the fact that mapping of 
fermionic spin glass transitions to Pseudo Yang Lee edge singularities have been proved may be enough motivation 
for an analysis of onedimensional cases. 

While it is complicated to solve ID fermionic Ising spin glasses ( no matter whether d = 1 is below the lower critical 
dimension or above due to long range interaction) exactly, the clean fermionic Ising chain offers simple exact solutions. 
These solutions discussed here will provide some insight into the role of the chemical potential and moreover generalize 
known resultsjato the complex /u-plane. 

Yang and LeeE3 derived the distribution of zeroes of the partition function of finite and infinite Ising chains within 
the complex magnetic field plane. In fermionic Ising systems the chemical potential can be seen as complemen- 
tary to the magnetic field. Stimulated by the representation of conventional Ising chains by fermionic ones with 
special imaginary chemical potential, one may wish to extend the Yang Lee analysis to a fourdimensional space 
((Ren), Im(n), Re(h), Im(h)). 

The transfer matrix T f of the fermionic Ising chain in a (complex) field h and with (complex) chemical potential /i 
reads 



O 0V 



1 



1 

-P/J, 



e ±0{»+h) e \0(h-p,) 



-0J 



-0J e l3(J-h) J 



(47) 



while the corresponding one of the S — ±l-chain is given by 



T s = 



-0J 



-0J p 0(J-h) 



(48) 



The transfer matrices and their eigenvalues don't map onto each other at /i = i%T, but the partition function does 
according to 



Z (N) = TtT n = (2i) N Z^ (49) 

for any number N of sites. 

The largest eigenvalue determines the free energy of the infinite chain, while the second largest is required in addition 
to determine the correlation length. The largest eigenvalue results in general from a cubic equation resulting in a 
lengthy result for the free energy. For several purposes it is sufficient to know the eigenvalues for (h = 0, /i), /x arbitrary 
complex. The half-filled case (fx = 0, h), h arbitrary complex, may also be considered separately. 
The eigenvalues for h = are found as 



A± = eP»[ch(0ii) + ch(0J) 

±yj(ch(l3n) + ch(f3J)) 2 + 4ch(/3fj.)(l - ch(f3J))} 
A = , Ai = 2e^sh(f3J). (50) 

The correlation length is given by 



17 



£ = l/j„(*±). (51) 
M 

In the T — > 0-limit a transition (formally) arises at /x = J and due to the properties 

Ai - ea;p(/3( A1 + J)) (52) 

/ exp{p(n + J)), fj,<J 
+ ~ \ exp(2pn), /u > J 

Thus the correlation length diverges only for /i < J, which is comprehensible since the energy required for adding a 
fermion is larger than the gain from a magnetic bond if /U > J. One obtains 

/ ezp(/3(2J - M )) ,0<^<J . . 

4 \ exp(/3J/2 , M = J. ^ 

while £ <~ T/(/i — J) for /i > J. Taking the derivative of the free energy / = —Tln\ + w.r.t. the chemical potential 
yields the filling factor showing that in the zero temperature limit the system is completely filled provided /i is positive 
(empty for negative n). Thus there is no physical T — O-transition of this simple system. The correlation length 
diverges in the T — > limit for all fillings v. 

The zero-field partition function shows that Yang-Lee zeroes approach = ±J for T — > 0. This means that for 
(h = 0, T = 0) nonanalytic behavior (as a function of the real chemical potential) can only occur at the values \x — ±J. 
It is instructive to consider N = 2 explicitly, which yields 

z n=2 = teWp^chtfrf + c h(/3J)f (54) 
+ch 2 (f3h){e 2f)J - 1) - sh(2/3J)} 
This almost trivial case already shows zeroes at 

fi = ±(J+{\ln2±i^)T) (55) 

while allowing for finite complex magnetic field the first zero different from ± J in the T — > 0-limit becomes possible 
with 

li = ±(i^ + 2imir)T = fj, + imwT. (56) 

More zeroes occur on the T — 0-axis as N is increased. In the N — > 0-limit a density function is expected in accor- 
dance with the divergent behavior of the correlation length for any \i. Its evaluation is beyond the scope of this paper. 



VII. MAPPING THE TWODIMENSIONAL ISING MODEL WITH IMAGINARY MAGNETIC FIELD 

H = IfT INTO FERMIONIC SPACE 

The complementary role of complex magnetic field and complex chemical potential can nicely be seen by recalling 
that the 2d Ising model with nonreal magnetic field had only been solved exactly at the special value H = i^nT. This 
value corresponds to /i = i^irT, which maps the fcrmionic Ising model onto the one above. Thus the exact solution 
of the 2d fermionic Ising model with /i = H = i\kT is known. Moreover this special model maps onto an interacting 
model which describes the interaction of spinless fermions with a special species obeying bare Bose statistics but 
carrying along the minus signs of fermion interactions. The Hamiltonian of this model can be written 

H = - J v^ a j - » S("*T + na) - h < (57) 

ij i i 

with n = h = i\T . While this value of fi provides the desired map between 2-state- and 4-state-per-site-space for 
any magnetic field, the special value considered here reduces to zero the imaginary field of one fermionic species, while 
the other acquires a shift equal to the distance between Bose- and Fermi-Matsubara energies. 
Setting c = a-f and d = aj, the Hamiltonian reads 
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cjcidjdj - d\didjCj + d\did^dj] 



iirT^d\di. 



(58) 



The effect of the new imaginary chemical potential iirT acting only on the d-fcrmions can both be imagined in 
standard many body diagram theory and within the fermionic path integral representation. In both formalisms the 
fermionic Matsubara energies of the d-fermions become bosonic ones. In other words d-propagators become bosonic 
while c-propagators remain (unchanged) fermionic. All interactions (vertices) between c-fermions and d-particles 
remain fermionic in character, ie they carry the minus-signs due to the anticommutation rules of the operators c 
and d. In the fermionic path integral representation Grassmann-fields are indispensable for both c- and d-operators. 
Denoting them by tPcIt) and by ipd(T~) respectively, the important change due to the z7rTn<j-term, which we absorb by 
the phase transformation exp(iirT)iJ)(T) — ip(r) in the new anticommuting fields tp, occurs in the unusual periodicity 
of tpdir) on the imaginary time axis: 



while iP c (t) = —iPc{t + P) retain the fermionic antipcriodicity. The periodic behavior of the ipd's purports of course 
into the bosonic propagator form. 

So far the mapping has been exact. In a perturbatively exact way one can even go beyond the above conclusion by 
stating that the 2D Ising model for this special imaginary magnetic field maps onto a coupled Fermi-Bose system 
(c-d) with the speciality that any d-loop contributes additional (— l)-factors. Despite the speciality of identical 
coupling constants between all species this system is physical in the sense that there are no imaginary fields, but 
the mapping of the known solution of the 2d Ising model discussed above is of order by order character. 
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